FINITE-SIZE SPIN- WAVE THEORY OF A COLLINEAR 
ANTIFERROMAGNET 



LUCA CAPRIOTTI 



Kavly Institute for Theoretical Physics, University of California at Santa Barbara 
Santa Barbara CA 93106-4030, United States of America 

caprio@kitp.ucsb.edu 



Received Day Month Year 
Revised Day Month Year 



The ground-state and low-energy properties of the two-dimensional J\ — Ji Heisenberg 
model in the collinear phase are investigated using finite-size spin-wave theory [Q. F. 
Zhong and S. Sorella, Europhys. Lett. 21, 629 (1993)], and Lanczos exact diagonaliza- 
tions. For spin one-half — where the effects of quantization are the strongest — the spin- 
wave expansion turns out to be quantitatively accurate for J2/J1 > 0.8. In this regime, 
both the magnetic structure factor and the spin susceptibility are very close to the spin- 
wave predictions. The spin-wave estimate of the order parameter in the collinear phase, 
tnt ~ 0.3, is in remarkable agreement with recent neutron scattering measurements on 
Li 2 VOSi0 4 . 
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1. Introduction 

The experimental realization of quasi-2d frustrated antiferromagnets, 1 such as 
Li2VOSi04, Li2VOGeC>4, and VOM0O4 has recently generated a renewed inter- 
est in the physics of the so-called collinear antiferromagnets. In these compounds 
the relevant superexchange interactions involve s = 1/2 V i+ ions on weakly cou- 
pled stacked planes, and the magnetic behavior is likely to be described by the 
J1 — J2 Heisenberg model on the square lattice, 



where Ji and J2 are the (positive) nearest- neighbor (n.n.) and next-nearest-neighbor 
(n.n.n.) couplings, respectively. In the experimentally relevant case J2 > J%, frustra- 
tion is known to originate in the classical (s — > 00) limit a collinear low-temperature 
phase. 2,3,4 In fact, for J2/J1 > 0.5, the ground state for s — > 00 is a state where the 
spins are ferromagnetically aligned along one direction and antiferromagnetically in 
the other, corresponding to magnetic wave vectors Q = (tt,0) or Q = (0, tt). These 
two families of states break not only the SU(2) symmetry and the translational in- 
variance of the Hamiltonian, as the conventional Neel state, but also its invariancc 
under n/2 real-space rotations. Interestingly, this additional two- fold degeneracy is 
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expected to generate non-trivial finite-temperature properties. In fact, an Ising-likc 
order parameter can be denned in order to discriminate between the two classes of 
states which are connected by real space ir/2 rotations, and a finite-temperature 
second-order phase transition is expected. 4 

How these low-temperature properties are affected by quantum fluctuations, es- 
pecially for the experimentally interesting s = 1/2 case, is an open problem currently 
preventing a faithful comparison with the experiments. 1 ' 5,6 In fact, due to numerical 
instabilities induced by frustration (sign-problem), 7 this issue cannot be investigated 
with the powerful stochastic numerical methods successfully employed for the n.n 
Heisenberg model 8 , so that approximate approaches have to be followed. 7 ' 5,6 

In this paper, we address the reliability of the spin-wave (sw) theory as an an- 
alytical tool to give quantitative predictions on the low-energy properties of the 
J\—Ji Heisenberg model in the collinear phase. To this purpose, a previously in- 
troduced finite-size sw theory 9 is generalized to the collinear phase and a direct 
comparison with Lanczos exact diagonalization results 10 is performed for s = 1/2. 
This allows us to demonstrate the effectiveness of sw theory in describing the low- 
energy properties of the J1 — J2 Heisenberg antiferromagnet in the experimentally 
relevant regime J2 > J\. 

2. Finite-size spin-wave theory 

The finite-size sw theory of Zhong and Sorella, 9 is a rigorous generalization of the 
standard sw approach to finite clusters with N sites, which allows one to avoid the 
spurious Goldstone-mode divergences - related to the SU(2) symmetry breaking 
assumption - in a straightforward way, without imposing any ad hoc holonomic 
constraint on the sublattice magnetization. 11 In the collinear case, a systematic 
finite-size sw expansion can be derived for J2/J1 > 0.5 considering as the approxi- 
mated ground state the symmetric linear combination 

hk«,> = (M> + |0,7r»/>/2, (2) 

where \n, 0} and |0, n) are the (normalized) sw ground states obtained assuming 
the classical collinear order with Q = (w, 0), and Q = (0, it), respectively. Since on 
any finite size the ground state is expected to have all the spatial symmetries of the 
Hamiltonian it is important to consider the symmetric linear combination (2) in 
order to obtain an accurate description of its correlations. In addition, this is also 
crucial to generalize rigorously the finite-size sw theory to the collinear phase. 

Both \ir, 0) and |0, n) in Eq. (2) can be obtained according to the standard sw 
procedure. The first step is to apply the unitary transformation which defines a 
spatially varying reference frame pointing along the local direction of the spins in 
the classical state, assumed in the xy plane: U = exp[7r/2^\(cos<3 ■ Ti — 1)5?]. 
Then, using the Holstein-Primakoff spin-boson transformation to first-order in l/s, 
S x = s — a^a, S v = ■\/|'(a t + a), S z = i^fj(a^ — a), and translational invariance, 
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the leading term of the Hamiltonian in the sw expansion reads 

H sw = E c i + 2 Jxs ^2 [A k ala k + -Bfc(a£al fc + a k a- k ) (3) 

k 

where E cl = —2J 2 Ns 2 is the classical ground-state energy, and a\ = 
l/y/N ^2 r exp[— ik ■ r]a\, with k = (k x , k y ) varying in the first Brillouin zone of the 
lattice . Here A k = 2 J 2 / Ji + cos k x and B k = — (cos k y + 2 J 2 / J\ cos k x cos k v ) for Q = 
(0, tt) and with k x *— > ky for Q — (tt , 0). Such leading part of the Hamiltonian, being 
bilinear in the Bose operators, can be diagonalized for k ^ (0, 0), (0, tt), (tt, 0), (tt, tt) 
using the well-known Bogoliubov transformation, a k = u k a k + v k c^_ k , with u k — 
^(A k + u k )/2u k , v k = -sgn(B k ) v / (A k -u k )/2u k , u> k = A\-B\ being the sw 
dispersion relation. In contrast, the Goldstone modes cannot be diagonalized by this 
transformation since u k and v k are not defined for the fc-vectors for which uj k = 0. 
For infinite systems such divergences are integrable in two dimensions and do not 
lead to any contribution. On any finite-size system, instead, they are important and 
they must be treated independently. 

To this end, it is possible to define a set of Hcrmitian operators commuting with 
each other and with the Hamiltonian: Q k = i(a k — a k ) for k = (0, 0) and (it, tt) and 
Q k = (a\ + a k ) for k = (0, tt) and (tt, 0). With these definitions the Goldstone- mode 
contribution to (3) reads 7i sm = Jis y~] s k m A k [Q^ — l] , with the sum extended 
only to the four singular modes. This allows us to diagonalize the Hamiltonian 
for any finite size in a basis where the Q k s have definite quantum numbers. In 
particular, being A k > for J2/J1 > 0.5, H sm favors a ground state with Q k = 
and it is easy to show that this implies to order 1/s S^OjO) = 5*^(0,0) = and 
S y (TT,Tr) = S z (tt,tt) = where S a (q) = J2 r ( S aSr)e iq ' r is the a component of the 
magnetic structure factor. These relations together with the sum rules ^(O, 0) = 
S x (tt,tt) = 0, following from the translational invariance of the sw Hamiltonian, 
are consistent with a singlet ground-state with a zero value of the total spin on 
each sublattice. This is rigorously true according to the Lieb-Mattis theorem 12 for 
J2/J1 — > 00. In fact, in this limit the system decouples into two independent n.n. 
Heisenberg models on the two sublattices. In general, a singlet ground state is also 
expected for the J1—J2 model on any finite size with an even number of sites and 
non-frustrating boundary conditions. The sum rules S a (n 7 Tr) = 0, instead, are not 
necessarily true on a finite size and for finite values of J2 / J\ but they turn out to 
be fulfilled with good approximation in the collinear phase (see below). 

According to this procedure, the sw estimate of the ground-state energy reads, 
E = E c i + Ji s J2 k [ w fe ~ Ak], and it turns out to be in remarkable agreement 
with the exact results for s = 1/2 and N = 36 in the whole range J2/J1 > 0.7 
(Fig. 1), i.e., far enough from the classical critical point J2/J1 — 0.5 separating the 
collinear from the Neel zero-temperature phases. In fact, in this strongly frustrated 
regime, a non-magnetic ground state of purely quantum mechanical nature is likely 
to be stabilized for s = 1/2 (Ref. 13). It is therefore natural to expect sw theory 
to fail in describing its correlations. Indeed, approaching the fully frustrated point 
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J2/J1 = 0.5, the sw expansion eventually breaks down leading to a zero value of the 
antiferromagnetic order parameter 2 and of the uniform susceptibility (see below) . 




A little more involved calculation allows us also to access the spin-rotation invari- 
ant spin-spin correlation functions: C(r) = (Si • S i+r ). The different contributions 
coming from the two degenerate states in Eq. (2), Cg(r), can be found by adding 
to the Hamiltonian a term of the form 7i{h r ) = J 2 /i r /2 J^. • Si +T where the 
vectors r are the z T equivalent lattice vectors that can be reached from a site i with 
|t| = \r\. The resulting sw Hamiltonian has the same form of Eq. (3) with field- 
dependent coefficients and it can be diagonalized similarly to the zero-field case. 
Then, according to the Hellmann-Feynman theorem, 14 the spin-spin correlations 
can be obtained by differentiating the sw energy in presence of the perturbation 
with respect to h r , after a careful treatment of the singular modes. In particular, in 
the collinear case, such rigorous treatment is possible only assuming for the ground 
state the symmetric form (2). 

The resulting correlation functions have different expressions depending on the 
form of the vector r. In particular, for r connecting the same sublattice, they read: 



C Q (r) = s 2 + sS rfi - — ^ (1 - cosfc • r) ^ , (4) 
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if e iQ r = 1, and 

1„(r) = -s 2 - 

N N 



a ns 

CqW - -s 2 - -T7 + —^2(v 2 k -cosk-ru k v k ) (5) 



k 



if e lQ r = —1, with the sums extended to all the non singular (ns) fc's. For r 
connecting different sublattices, instead, the spin-spin correlations are: 

ns 

Cq{r) = — J2( T * v k ~ f kU k v k ) , (6) 
k 

where = l/z T ^2 T cos A; • r with e l( ^' T — 1, and t\ = l/z T J2 T C0S ^" T 
with e l< 5' T = — f. The sw spin-spin correlation functions calculated as C(r) = 
(C(o,7r)( r ) + C(7r.o)( r )) /2 are invariant under the symmetry operations of the lat- 
tice, as expected for the ground-state expectation values on any finite size. Using 
these expressions the spin-isotropic magnetic structure factor, S(q) = J2 r C{r)e lq ' r , 
can be also calculated. In addition, having obtained an ordered expansion in 1/s 
for the spin-spin correlation functions, C(r) = s 2 ( e 4 ( 7r <°)- r + e 4 (°' 7r )' r )/2 + sa(r), 
the antiferromagnetic order parameter, m) = y / 2/NS(Q), 15 can be expanded as 
m) = s(l + a/s) with a = 1/Nj2 r e lQ ' r a(r), and it reads: 



ns 

1 \ ^ 2 

Ns Ns 4^ Vk 



(7) 



Note that for — > oo the known sw result 2,3 for the thermodynamic order param- 
eter, m) = \J (0, 7r|(5fp|0, 7r), is recovered. As shown in Fig. 2, the sw magnetic 
structure factor, S(q), is in remarkable quantitative agreement with the exact di- 
agonalization result for s = 1/2 and J2/J1 = 1 for all the wave vectors q, and the 
sw sum rule S(it 7 it) — is fulfilled almost exactly. In particular, the sw approach 
provides a very accurate estimate (within a few percent) of the antiferromagnetic 
order parameter in the whole range of J2/J1 > 0.7. 

The sharp drop of the order parameter at J2/J1 — 0.7 is suggestive of a first 
order transition between the gapped and the collinear phase, as also indicated by 
a recent variational study. 16 This leads in particular to a very fast saturation of 
the order parameter to the value of the n.n. Heisenberg model, 17 expected for 
J2/J1 —> 00. Remarkably, the sw prediction for J2/J1 > 1.0 is in excellent agreement 
with the result of recent neutron scattering measurements 18 on Li2VOSi04, giving 
m) = 0.31(2) in the zero temperature limit. This compound is known to have 
Ji ^ -Ai even if the precise value of the frustration ratio is presently much debated. 
1,6,19 Unfortunately, the weak dependence of the order parameter on J2/J1 for 
J2/J1 ^ 1-0 does not allow to use the sw prediction to determine more precisely 
this ratio. 

Within the finite-size sw theory, it is also possible to study the structure of 
the low-energy spin excitation spectrum which is connected to the ground-state 
antiferromagnetic correlations and to the uniform susceptibility. Following Lavalle 
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Fig. 2. Lower panel: magnetic structure factor for J2/J1 = 1, and N = 36. Spin-wave (full dots 
and dashed line), exact (empty circles and dotted line). Upper panel: antiferromagnetic order 
parameter at wave vector Q as a function of the frustration ratio J2/J1, for N = 36. Spin- wave 
(dashed line), exact (empty circles). The arrows and the full dot indicate the sw and the exact 
17 results for the n.n. Heisenberg model, respectively, and the star is the experimental result for 
Li2VOSi04 18 which has a still undetermined J2/J1 > 1. 1 ' 6 > 19 The continuous lines in both panels 
are the spin- wave results in the thermodynam ic limit. 



et al, in order to stabilize an excitation of total spin S, a Zecman term, 7i(h) = 
~hs^2 i Sf , is added to the spin Hamiltonian. Classically, for magnetic fields not 
large enough to induce a spin-flop transition, the classical solution is simply obtained 
by canting the spins of an angle Oh along the direction of the field h [sin#/j = 
/i/4Ji(l + 2J2/J1)]. After a new rotation of the reference frame in order to align 
it to the new directions of the spins, the finite-size sw expansion is straightforward 
and again it leads to a linearized Hamiltonian of the same form of Eq. (3) with field- 
dependent coefficients A\ = Ak + Bk sin 2 9h, B% = £?fc(l — sin 2 Oh), and classical 
energy. In this case, due to the reduced spin-rotation symmetry of the Hamiltonian 
only the (0, 0) and the (it, it) modes are singular. The former favors a value of 
S z (0,0) = N(ss'm6) 2 consistent with the applied field, at the classical level, while 
the latter generates the usual approximated sum rule S z (ir, it) — 0, which is not 
affected by the uniform field along the z direction. The total spin S = N(Sf) of the 
excitation can be related to the magnetic field h by means of the Hcllmann-Feynman 
theorem 

1 d h 
{S > )h = -N- S dh m = U M l + 2J 2 /J 1 ) X 
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2iY S (l + 2J 2 M) V l A k~ B k 



(8) 



where 



E(h) = E cl - (shf ^^^ + K - 4] , (9) 



and uu k = y (A^) 2 — (£>£) 2 . Using these expressions the energy spectrum E(S) can 
be calculated by means of a Legcndrc transform, E(S) = E{h) + hsS. Besides, the 
uniform susceptibility in the thermodynamic limit, x = — l/Ns 2 d 2 E(h)/dh 2 \h=o, 
can be obtained from Eq. (9) by direct derivation and it reads: 



/ 1 1 f d 2 fc R / A k + B k 

XlXd = 1 " 2 S (1 + 2J 2 M) J \ ~Ak~-B~k • (10) 

where Xcl = l/4Ji(l + 2J 2 / J\) is the susceptibility in the classical limit. 

Due to the discreteness of the energy spectrum, a direct comparison of the sw 
susceptibility with the exact diagonalization results is not possible on a finite size. 
However, in a quantum antiferromagnet with a long-range ordered ground state, 
the susceptibility is directly related to the properties of the low-energy spectrum 
and this can be exploited to establish the accuracy of the sw predictions. In fact, 
whenever long-range order is present in the thermodynamic limit, the low- lying 
spin excitations (S -C \/~N) order themselves as in the spectrum of a free quantum 
rotator, E(S)—E — S(S+1)/2IN, where I is the so-called momentum of inertiaper 
site. This definition is very close to the definition of the uniform spin susceptibility, 
which can be calculated by taking first the infinite-volume limit of the energy per 
site at fixed magnetization m = S/N, e(m) = E(S)/N, and then letting m — > 0: 
e(m) = e + m 2 /(2x). As it is known from the low-energy theory of a quantum 
antiferromagnet, 21 this allows for an identification between / and \/N so that 
the quantity [2%sP = N (E(S) — E ) [S(S + 1)]~ approaches for infinite size the 
physical inverse susceptibility 1/2%, for any spin-excitation S -C N. This features 
are clearly verified by the sw scheme as it is illustrated in the lower panel of Fig. 3 
for J2/J1 = 1, where l/2%s is plotted for S = L = ^/N and approaches the 
thermodynamic value calculated according to Eq. (10) (l/2x sw — 11-36). 

The comparison with the exact results of the sw spectrum E(S) — Eq vs S(S+ 1) 
(inset of Fig. 3), and, equivalently, of the quantity l/2xs=i (upper panel of Fig. 3) 
demonstrates the accuracy of sw theory in estimating the spin susceptibility of the 
J1 — J2 model in the collinear phase for Jij J\ > 0.8. As it is shown in the same 
figure (upper panel), the classical uniform susceptibility is strongly renormalizcd 
by quantum fluctuations for s = 1/2 at the sw level. As expected, approaching 
J2/J1 = 0.5 such reduction is enhanced due to the increased frustration and it 
leads eventually to the vanishing of the susceptibility for J2/J1 — 0.507 and to 
the expected transition to a spin-gapped non-magnetic phase. 7 However, from the 
comparison with the exact diagonalization results for both the uniform susceptibility 
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and the order parameter, it appears clear that the effects of quantum fluctuations 
are underestimated by the sw expansion in the regime of strong frustration so that 
the transition to the non-magnetic region is likely to occur at slightly higher values 
of the J2/J1 ratio (J2/J1 — 0.6). 




0.1 0.2 0.3 

1/L 

Fig. 3. Lower panel: size dependence of the spin-wave 1/2\l (see text) for J2/J1 = 1. Inset: spin- 
wave (full dots and continuous line) and exact (empty circles and dotted line) low-energy spectra 
as a function of |S 2 | = S(S + 1) for N = 36, and J2/J1 = 1. Upper panel: uniform susceptibilities 
as a function of J2/ 'J\. Dashed line: spin- wave XS=l f° r N = 36. Empty circles: exact Xs—l f° r 
N = 36. Continuous line: thermodynamic spin- wave \. 



3. Conclusions 

In conclusion, we have presented a finite-size spin- wave study of the J1 — J2 Heisen- 
berg antiferromagnet in the collinear phase, focusing in particular on the ground- 
state spin-spin correlations and the low-energy excitations. For s = 1/2, the 
comparison with exact diagonalization results reveals a remarkable agreement for 
J2/J1 > 0.8. In particular, both the antiferromagnetic structure factor and the spin 
susceptibility are very close to the spin-wave predictions. The accuracy of the re- 
sults on finite sizes indicates that the spin-wave expansion provides a quantitatively 
reliable description of the ground state of the J1—J2 model in the experimentally rel- 
evant regime J2 > Ji . This is confirmed by recent neutron scattering measurements 
on Li2VOSi04, providing a value of the order parameter in excellent agreement with 
the spin-wave prediction in the collinear phase, ~ 0.3. 
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The accuracy of spin-wave theory in describing the low-energy excitations also 
suggests that theoretical approaches to the thermodynamics based on a Gaussian 
treatment of quantum fluctuations, such as effective Hamiltonian methods, are 
likely to successfully describe the finite-temperature properties of the Ji — J2 model 
in the collinear phase. Work on this line of research would allow to establish a 
closer contact between the frustrated Hciscnberg model and the physical proper- 
ties of the recently synthesized collinear antiferromagnets, helping in their precise 
characterization and in the determination of the actual frustration ratios. 
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